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NAME:  Dr.  Arje  Nachman 

BRIEF  DESCRIPTION  OF  PORTFOLIO: 

Interrogation  (Modeling/Simulation)  of  Linear/Nonlinear 
Maxwell’s  Equations 

LIST  SUB-AREAS  IN  PORTFOLIO: 

Theoretical  Nonlinear  Optics 
Wave  Propagation  Through  Complex  Media 
Fundamentals  of  Antenna  Design/Operation 
Fundamentals  of  Effects  of  EM  Exposure  on  Circuitry 
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Scientific  Challenges 


•  Wave  Propagation  Through  Complex  Media 


Details  of  time-domain  dynamics  of  EM  pulses  through  Dispersive, 
Conductive,  and/or  Random/Turbulent  media 

Research  provides  optimism  that  the  class  of  waveforms  termed  Precursors 
have  the  potential  to  upgrade  imaging  quality. 

•  Antenna  Design/Operation 

Suitable  PARTNERSHIPS  of  MATERIALS  and  GEOMETRY  can  deliver 
man-made  composites  which  exhibit  novel  EM  attributes. 

Such  METAMATERIALS  include:  NIMs,  PBGs,  and  “Unidirectional” 
composites. 

Growing  reliance  on  small  UAVs  drives  the  need  to  miniaturize  antennas  and 
make  them  more  responsive. 
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Scientific  Challenges 


•  Nonlinear  Optics 

Fundamental  modeling/simulation  research  which  addresses  concerns  with 
femtosecond  filament  arrangements  and  plasma  channel  characteristics. 

Advances  in  modeling/simulation  of  fiber  and  solid  state  lasers  to  guide  the 
development  of  compact,  high  energy  systems. 

•  RF  Effects  on  Circuitry 

Identification  of  waveforms  which  produce  various  realizations  of  circuit  upset 
(includes  chaos). 

Complicated  by  the  fact  that  effects  are  likely  to  be  dictated  by  the  activity  of 
the  circuit  (eg,  routines  being  run  by  laptop). 
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MURIs 


This  portfolio  has  an  existing  MURI  (Ultrashort  Laser  Pulses) 
which  just  completed  its  1st  year. 

Two  more  MURIs  supportive  of  the  portfolio  subarea 
“Wave  Propagation  Through  Complex  Media” 

will  be  starting: 

•  “Deep  Optical  Turbulence  Physics” 

•  “High  Power,  Low-Loss,  Artificial  Materials  for 
Transformational  Electromagnetics” 
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\  A  Mathematical  Modeling  and  Experimental  Validation 
of  Ultrafast  Nonlinear  Light-Matter  Coupling  Associated 
with  Filamentation  in  Transparent  Media 

MURI  PMs— Dr  Arje  Nachman  and  Dr  Enrique  Parra 
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New  State  of  Matter:  Uncorrelated  Electrons 

to  Plasma  Transition 

B.  Pasenow,  M.  Brio,  J.V.  Moloney,  S.  W.  Koch,  S.H.  Chen,  A.  Becker,  and  A.  Jaron-Becker 

•  MURI  collaboration  between  University  of  Arizona  and  University  of  Colorado 

•  USP  ionization:  highly  anisotropic  uncorrelated  electron  distributions 

•  Relaxation/isotropization  of  non-equilibrium  distributions  due  to  Coulomb  scattering 

•  THz  EM  response:  evolution  of  loss  of  anisotropy  and  carrier  density  (plasmon  pole) 
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Fermi-Dirac  electron  momentum 
distribution  in  isotropic  plasma 
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^  FORMATION  of  PRECURSORS 

and  their  ALGEBRAIC  DECAY  with  DISTANCE 

into  DISPERSIVE  MEDIA 


Dr  Kurt  Oughstun/UVt  and  Dr  Natalie  Cartwright  (YIPVSUNY-NewPaltz 


AIR 


Electric  field  component 
before  incidence: 

SQUARE-WAVE  MODULATED  SINE 


Note:  The  CW  version  of 
the  above  sinusoid  would 
experience  exponential 
decay!  Z  =  0 


MATERIAL 

water,  etc 


HUH 


Electric  field 
component  of  the 
transmitted  pulse. 


Electric  field  component 
of  the  propagated  pulse. 
This  precursor  decays 
as  the  inverse  square 
root  of  z  and  contains  a 
generous  bandwidth. 


Z>0 
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THE  PROPAGATED  PULSE 

how  they  did  it 


The  electric  field  component  of  the  propagated  pulse  on  any  plane  z  >  0  is  given  by 


1  f  ia +oo 

E(z,t)  =  —  Rel  i  J  Te(cd)E(z  <  0,<o)e 


ik{co)z-icot 


In 


ia—  oo 


dco\. 


~  co 

where  k(co)  =  —  /z2  (&>)  is  the  complex  wave  number  of  the  dielectric  material. 

TE  O)  = ; — n2  O)  = 


1  +  772(&>) 


1 


l- 


-atn+liSco 


For  values  of  ct/z=  1  the  contour  may  be  closed  in  the  upper  half  plane  and  application  of 
Jordan’s  lemma  gives  /  )  =  O,  /  <c  z  /  c. 

This  integral  representation  of  the  field  has  no  known  exact  solution  when  ct/z>1 


Asymptotic  methods,  such  as  saddle  point  methods,  may  be  used  to 
find  an  approximation  to  the  propagated  pulse. 

These  methods  require  the  deformation  of  the  Bromwich  contour  through  the 
valleys  of  the  accessible  saddle  points  of  the  complex  phase  function 

(j){co, 6)  -  ico[n(co)  -  0\  0-  — 


which  is  completely  characterized  by  the  dielectric  material. 
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ACTIVE  INFRARED  IMAGING  THROUGH 
SPARSE  DISCRETE  MEDIA 


Drs.  Elizabeth  and  Marek  Bleszynski 
Monopole,  Inc 


PROJECT  OBJECTIVE 


investigate  possibility  of  using  wide-band,  infrared  (IR)  pulses 
in  imaging  through  media  composed  of  sparsely  distributed 
discrete  particles ,  such  as  clouds ,  fog,  dust,  or  smoke 


“ sparse ”  medium:  mean-free  path*  much  larger 
than  average  separation  between  scatterers 

*mean  free  path  is  the  penetration  depth  over  which  intensity  decreases  by  1/e 
so  sparse=mean  free  path  corresponding  to  the  highest  frequencies  in  the  pulse 
spectrum  is  much  larger  than  average  separation  between  scatterers. 
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PROPAGATION  OF  A  WIDE-BAND IR  PULSE 
THROUGH  SPARSE  DISCRETE  MEDIA 


example:  a  trapezoidally  modulated  IR  pulse  propagating  through  a  cloud 


pulse  parameters:  cloud  parameters: 

-  carrier  wavelength:  A  =  10pm  (v  =  30THz)  -  average  droplet  radius:  a  ~  5|jm 

-  number  of  cycles  in  the  pulse:  m  =  60  -  average  droplet-droplet  distance:  R  ~  1  mm 

-  rise/fall  time:  0.05  of  the  carrier  period 


transmitted  pulse 
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i PROPAGATION  OF  A  WIDE-BAND IR  PULSE 
THROUGH  SPARSE  DISCRETE  MEDIA 


-  medium  as  a  filter  attenuating  high  frequencies 

-  evolution  of  the  propagating  pulse  into  a  Brillouin  precursor 

pulse  spectrum  after  propagating 


pulse  amplitude  after  propagating 
z  =  100m,  200m,  300m 
through  the  cloud  medium 


z  =  100m,  200m,  300m 
throuah  the  cloud  medium 

— — . . . • — • — • — i — 


> 


Brillouin  precursor-type  structures 
associated  with  the  leading  and 
trailing  edges  of  the  transmitted  pulse 


100  1000  10000  100000  100000C 

1  THz  V  (GHz) 

the  transmitted  pulse  spectrum  has  a  significant 
amount  of  low-frequency  components 
the  high-frequency  part  of  the  spectrum  is 
attenuated  very  strongly  as  the  pulse  propagates 
the  low-frequency  part  of  the  spectrum  isl ^ 


weakly  attenuated 
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! PROPAGATION  OF  A  “CHIRPED  TRAIN  OF  PULSES 


generate  a  coherent  train  of  N  pulses  (with  small  rise/fall  times) 
emitted  at  linearly  increasing  pulse-pulse  intervals 


linearly  increasing  time  intervals  T  between  pulses 

r 

chirped  train  characterized  by:  effective  center  frequency  and  effective  bandwidth 

1  1 


1 


= 


(r) 


B  c  = 


T 


max 


after  passing  through  clouds  the  pulse  train  becomes  a  train  of  precursor-type  pulses 
associated  with  leading  and  trailing  edges  of  the  transmitted  trapezoidal  pulses 
and  attenuated  approximately  algebraically  (not  exponentially) 


t 
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PRECURSOR  SUMMARY 

<*» 

•  For  MW  (radar)  imaging,  it  is  difficult  to  produce  square-wave 
modulated  sinusoids. 

But  it  is  possible  for  modern  radars  to  produce  PRECURSORS! 

1  .These  interesting  waveforms  contain  greater  bandwidth  than 
conventional  narrow-band  radar  pulses, 

2.  Decay  algebraically  with  depth, 

3.  Experience  reduced  “flash”  at  any  media  interface, 

4.  Allow  for  easier  Matched  Filtering. 

•  None  of  this  is  seriously  modified  by  oblique  incidences. 

•  For  laser  imaging  through  clouds  (ladar)  it  is  very  easy  to 
produce  (nearly)  square-wave  modulated  sinusoids,  but  not 
precursors.  The  cloud  produces  the  precursor,  which  decays 
algebraically  with  depth.  In  order  to  compensate  for  the  loss  of 
the  high  IR  frequencies  (and  thus  detailed  spatial  resolution)  the 

notion  of  chirping  the  pulse  train  was  invented. 
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Laser  Beams  &  Turbulence 

Dr  Greg  Gbur/ElectricalEngineering 

Optical  beam  propagation  in  the  atmosphere  is  hindered  by  atmospheric 
turbulence:  random  fluctuations  of  the  refractive  index. 


Applications  such  as  free-space  optical  communications  and  LIDAR  are 
adversely  affected  by  atmospheric  turbulence,  which  induces  intensity 
fluctuations  (scintillations). 


Laser  beam  distorted  by  turbulence  ‘Heat  shimmer’  (mirage)  on  a  hot  day 
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^  ^Partially  Coherent  Beams  in 

'»* _ Turbulence _ 


>A  coherent  laser  essentially  propagates  its  energy  through  a 
single  coherent  beam,  which  may  partly  or  wholly  miss  detector 

>  The  partially  coherent  beam  sends  its  energy  through  multiple 
beamlets,  increasing  the  likelihood  of  hitting  a  detector  and 
smoothing  out  fluctuations 


>  Beamlets  need  to  be  significantly  different,  or  diverse 
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Airy  Beams  in  Turbulence 


A  spatial  filtering  system,  with  an  appropriate  phase  mask,  can  produce  an  Airy  beam 


source  plane  detector  plane 


Try  four  beamlets  designed  to  curve  together  at  target: 
each  beamlet  sees  different  turbulence  but  average  peak 

intensity  is  at  center  of  detector 
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On  axis  scintillation  index 


Scintillation  reaches  its 
minimum  possible  value 
with  only  four  beamlets! 

Evidently  four  beamlets  is 
“good  enough”  for  reducing 
intensity  variance! 

Additional  beamlets  are  not 
diverse  enough  to  provide 
additional  improvement 

Nonuniformly  polarized 
beams  another  idea 
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Imaging  with  Spatial  Nonlinearity 

Dr  Jason  Fleischer 
Electrical  Engineering,  Princeton 


•  Scalar  field  qj  obeys  nonlinear  Schrodinger  equation. 

•  Numerical  solution  via  split-step  Fourier  method. 


dyr 

dz 


i±-V2±+iAn(\y,\2) 

Linear  Nonlinear 


y/  =  \D  +  7V(I  y/  I  ) \y/ 


i//(x,y,z  +  dz ) 


=  ebdz+Adzy/(x,  y,  z) 
«  eDdzeNdzy/(x,  y,  z) 


Back-propagation  gives  yj(z)  if  we  know  qj(z  +  dz): 

y/(x,  y,  z)  ~  e~Ndze~Ddzy/(x,  y,  z  +  dz) 


Initial  value  problem 


Reconstruction 

Dimntinn 


Input 


Output 

(Measured) 
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Improved  imaging  with  NL: 
field  of  view  and  resolution 


Screen  blocks  rays. 


200  pm 


Sr0  75Ba0  25Nb206  Linear  Output 
Photorefractive  crystal 


Nonlinear  element 
retrieves  lost  bar 

30%  increase  in  FOV 
An  /  n0  ~  6  X 1 0-4 


Nonlinear  Output 


I  ft  I 

t  lit 


V linear-  0.095  VNL  -  0.32 


Airy 
Disks 


ffi?h 


Rayleigh  Criterion 

Visibility  Raylelgh  =  0.15 
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Nonlinear  Beam  Deflection  in 
Photonic  Lattices  with  Defects 


Optically-induced  photonic  lattice  with  defects: 
lattice  spacing  ~  10  microns 
defect 


1 D  lattice 


2D  lattice 


defect 


A  defect  is  a  local 
abnormality  inside  a  periodic 
lattice.  Experimentally 
created  1 D  and  2D  lattices 
with  single-site  defects  are 
shown  in  the  left  figures. 

Experimental  procedure  to  optically 
“write”  lattices  with  defects: 

(1)  let  a  laser  beam  pass  through  an 
amplitude  mask 

(2)  use  frequency  filtering  to  remove 
half  of  the  spatial  frequencies 

(3)  slightly  tilt  the  resulting  beam 


Dr  Jianke  Yang  (Math,  Univ  Vermont) 
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When  a  nonlinear  beam  is  launched  into  a  lattice  with  a  defect 
one  finds,  both  theoretically  and  experimentally,  that  at 
small  incident  angles  the  beam  is  reflected  by  the  defect  but  at 
large  incident  angles,  the  beam  passes  the  defect. 


small  angle: 
reflection 


large  angle: 
transmission 


z;  propagation  direction 

black  arrow:  probe-beam  direction  before  reaching  defect 
green  arrow:  probe-beam  direction  after  hitting  the  defect 


Result  shows  a  way  to  use  a  beam’s  incident  angle  to  control  its 
propagation  direction  in  a  lattice  network. 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


22 


Theoretical  Modeling 


.dU  1  ,d2U  d2U , 

i +  —  ( — r  + — T) 

dz  2kx  dx 2  cty2 


1,  3 

—  kji„  n 


W'e'33 


l  +  I(x,y)+\U\- 


■U  =0, 


/  =  /0cos2  — r(l  — g 


■x2  Id2 


d 


), 


ID  case 


/  =  /()  cos 


2  nx  zTiy 


d 


cos 


(1-e 


-( x2+y2)/d 2 


), 


2D  case 


U:  electric  field;  E0  :  applied  dc  field;  r33:  electro-optic  coefficient; 


70  :  lattice  intensity;  d:  spacing;  k{),k{ :  wave  numbers;  :  refractive  index 


Simulation  result: 

lattice  with  defect 


reflection  at 
small  angle 


transmission  at 
large  angle 
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Can  a  non-Hermitian  operator 
exhibit  real  spectra? 

_  volume  so.  number  24  physical  review  letters  isjunei99$ 


Real  Spectra  in  Xon-Hermitian  Hamiltonians  Having  T'T  Symmetry 

Carl  M.  Bender1  and  Stefan  Boettcher1' 

1  Department  of  Physics.  I Vashtngton  Untvmtiy.  St  Louts .  Missouri  6S130 
2  Center  for  Nonlinear  Studies,  Los  Alamos  National  Laboratory.  Los  Alamos,  .Yew  Mexico  S7S45 
’  CTSPS  Clark  Atlanta  Lnnerxth  Atlanta  Georgia  kOA14 

(Received  1  December  1997,  revised  manuscript  received  9  Apnl  1998) 

The  answer  is  yes  as  long  as  the  Hamiltonian  respects  PT^symmetry! 


P  \  x  — ^  — x  j  T  \  t  — ^  — t 

d *¥  fi  d2x¥ 

ifi  —  + - -  -  V (a)1?  =  0  PT  symmetry 

dt  2m  dx 


V(x)=V*(-x) 


Quantum  mechanical  oscillator 


PT  pseudo-Hermitian  oscillator 


T't symmetry  in 
Quantum  Mechanics  and  Optics 


W 


Schrodinger  Equation 


n2  d2v  ...  . w  . 

iti  —  + - —  -V(jc)T  =  0 


iX 


dt  2m  dx 

Time 


Paraxial  Equation 

dE  X2  d2E 

—  +  n(x)E  =  0 


+ 


dz  2  dx' 


t 


Planck’s  constant 


fi 


Propagation  distance  z 
wavelength 


X  =  - 
k 


Probability  amplitude  *¥(x,t) 


Electric  field  envelope  E(x ,  z) 


Complex  Potential  V ( x )  =  (x)  +  iVl  (x)  Complex  refraction  n(x)  —  nR  (x)  +  in1  (x) 

The  imaginary  part  ft7(x)  corresponds 

to  gain  (if  nT  <  0  )  or  loss  ( if  nT  >  0  ) 


So  nR  is  even 
and  n,  is  odd 


V(x)  = 


*  TT symmetry  condition 

V(-x)  n(x)  =  n(-x) 


* 
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w  ‘P'T symmetric  structures  and  devices 


PT  symmetry  in  optics  can  be  readily  established  by 
deliberately  involving  gain  and  loss. 


By  doing  so,  new  PT-symmetric  structures  and  materials  with 
useful  functionalities  can  be  envisioned. 


r»(Cr)=4. 272+4  92Si 


Loss-induced  transparency 
due  to  PT-symmetry  breaking 


On-chip  optical  isolators  and  circulators 


Dr  Tsampikos  Kottos  (Physics/Wesleyan)  and  Dr  Demetrios  Christodoulides  (Optics/UCF) 
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Electrically  Small 
Supergain  Arrays 


0.15  X 


0.15  X 


0.05  X 


Patent  Awarded -Dec  2011 

(Yaghjian,  Altshuler,  O’Donnell,  Best) 


*  Electrically  Small  (ES)  arrays  attain  7  dB  realizable  gain, 
5  dB  higher  than  any  previous  ES  antenna. 

*  Can  be  used  to  replace  much  larger  Yagi  antennas. 

*  Paves  the  way  for  3  or  more  element  endfire  arrays  with 
higher  gains  &  for  multiband  supergain  arrays. 


Future  research  to  increase  bandwidth  is  desirable. 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


27 


RF  Metamaterials  for  FOPEN  Application 

An  Emulation  9f  Anisotropy  Degenerate. Band  Edge 
Tx  une  Using  Standard  Printed  Circuit 


Lockheed  Martin  and  OSU 


Meta  FOPEN  Antenna  achieved  the  following  performance  (230  to  900  MHz): 

•  Bandwidth  >  3X  of  State  of  the  Art  (SOA) 

•  Half  size  of  SOA  (24"  x  24"  area  x  6"  height) 

•  Light  weight  (5.6  Lbs) 

•  Dual  polarized  with  <  20  dB  cross-pol  isolation 

•  High  power  capability  and  electronic  beam  steering 

•  Simple  and  low  cost  construction,  stamdard  PC  and  manufacturing  method 


Metamaterial  Equivalent  Circuit 
Degenerate  Band  Edge  (DBE) 


L-  Element  Conductor 
C- Overlapped  Element 
Shunt  L-  GP  Affect; 

Z  =  jZ0tan(kd) 

C-Element  above  GP 


Predicted  and  measured 
Beam  Steering 


Measured  Radiation  Patterns 


Scanning  Angle  20*  from  200  to  900MHz 
RnNzwJ  Gain  Pattern  PMOtteg 


integration  to  Shadow 
Harvest  Sensor  Pod 
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ELECTROMAGNETICS 
LAB  TASKS 


Dr.  Brad  Kramer(AFRL/RY),  “Electromagnetic  Materials  and  Antennas”  * 

1.  Model  Electromagnetically  Small  Antennas:  superdirective,  wide-band,  conformal 
Dr.  Ilya  Vitebskiy  (AFRL/RY)  “Metamaterials  for  the  Enhancement  of  Light-Matter  Interaction”  * 

1.  Performance  enhancement  of  various  transceivers 
Dr.  Saba  Mudaliar  (AFRL/RY),  “EM  Scattering  Studies” 

1.  Predict  scattering  from  clutter  and  rough  surfaces 
Dr.  Kris  Kim  (AFRL/RY),  “Predict  Far-Field  RCS  via  Near-Field  Data” 

(Dr.)  Jason  Parker  (AFRL/RY),  “Moving  Target  Radar  Feature  Extraction” 

Dr.  Nicholas  Usechak  (AFRL/RY),  “Dynamics  of  Reconfigurable/Agile  Quantum  Dot  Lasers”  ** 

1.  Investigate  control  of  amplitude-phase  coupling  in  Quantum  Dot  laser  systems 
Dr.  Timothy  Clarke  (AFRL/RD),  “Modeling  of  HPM  Effects  on  Digital  Electronics”  ** 

1.  Derive  mathematical  model  predicting  effects  (upset)  on  digital  electronics  when  exposed  to  various 
incident  EM  pulses 

Dr.  Danhong  Huang  (AFRL/RV),  “Models  for  Ultrafast  Carrier  Scattering  in  Semiconductors” 

1.  Model  IR  amplifier  for  extremely  weak  signals  and  distant  targets 
Dr.  Analee  Miranda  (AFRL/RY),  “Detection  and  Imaging  of  Underground  Facilities  Using  SAR  Data”  * 

Dr.  Matthew  Grupen  (AFRL/RY),  “Electronic  Band  Structure  for  High  Speed  Quantum  Electron  Device  Simulation" 
1.  Modeling/Simulation  of  quantum  tunneling  devices 

Dr.  lyad  Dajani  (AFRL/RD),  “Time  Dynamics  of  Stimulated  Brillouin  Scattering  in  Fiber  Amplifiers  with  Frequency 
Modulation” 

1.  SBS  suppression  research  to  realize  higher  power  in  narrow  linewidth  fiber  amps 
Dr.  Erik  Bochove  (AFRL/RD)  “Modeling  of  Large  Nonlinear  Passively  Phased  Fiber  Laser  Arrays”  ** 

*=New  for  FY12  **=Renewal  for  FY12 
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Subareas  Funding  Trends 


Wave  Propagation  Through  Complex*  Media 

*  Dispersive,  Conductive,  Random/Turbulent,  Man-Made  Composites 


Antenna  Design/Operation 


Effects  of  EM  Exposure  on  Circuitry 
Nonlinear  Optics 

(MURI  on  “Propagation  of  Ultrashort  Laser  Pulses  through  Transparent  Media”  began  1  Oct  2010) 


m 
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Connections  with 
Other  Organizations 


Extensive  interaction  with  Dr  Richard  Hammond/ARO  on 
UltraShort  Laser  Pulse  propagation  through  air 

Dr  Hammond  served  on  my  FY10  USLP  MURI  evaluation  panel  and 
I  served  on  his  FY11  USLP  MURI  panel 


•  JTO 

I  manage  the  JTO  MRI  “High  Power  Lasers  Using  Optically  Pumped 
Semiconductor  Laser  (OPSL)  Concepts”  which  ends  in  Aug  2012 

•  NRO 

•  Extensive  discussions/visits  regarding  impact  of  6.1  research  on  NRO  needs 

•  Arranged  for  2  Pis  to  participate  in  the  FY12  NRO  Seminar  series 
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Connections  with 
Other  Organizations 


•  ONR 


MURI  (U  Maryland)  “Exploiting  Nonlinear  Dynamics  for  Novel 
Sensor  Networks”  managed  by  Dr.  Michael  Shlesinger,  ONR 

I  serve  on  this  ONR  MURI  panel 

Negative  Index  Media  MURI 

Attended  review  of  ONR  (Dr.  Mark  Spector)  NIM  Metamaterials  MURI 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


